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Abstract 

We obtain and study an analytical solution of a de-Sitter thick domain wall in 
five-dimensional Einstein gravity interacting with a scalar field. The scalar field 
potential is axion-like, V{(f)) = a + 6cos(y^2730) with constants a,b satisfying 
—36 < 5a < 36, and the solution is expressed in terms of elliptic functions. 
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Solitons play major roles in non-perturbative aspects of field theory and string theory. 
Since string theory contains gravity, solitons necessarily couple with gravity in string theory. 
Because of the non-linearity and instability of the gravitational interactions, the inclusion 
of the gravitational dynamics into a soliton dynamics is a highly non-trivial problem. The 
simplest soliton interacting with gravity would be a domain wall in a coupled system of a scalar 
field and gravity. The easiest to treat are the BPS domain walls in supergravity theories, in 
which the equations of motions are given by the BPS first order differential equations [0. These 
solutions are static. In fact, it turns out that static domain walls without supersymmetry can 
also be treated in the same manner [@, ^ A solution of a domain wall with a de-Sitter 

expansion was first constructed in the thin wall approximation 0. Not so many analytic 
solutions of a thick domain wall with a de-Sitter expansion are known 0, In this paper, 
we obtain another analytic solution of a thick domain wall with a de-Sitter expansion. 

The system we discuss in this paper is the coupled system of gravity and a scalar field in 
five-dimensions, and the action is given by 

S = J dtd?xdy^ (^R - ]^g^^-d,<pd,cP - . (1) 

The metric ansatz we use is the warped geometry 

ds^ = a{y)\-de + e^^'Xdx'f) + dy\ (2) 

where H denotes the Hubble constant of the four-dimensional de-Sitter space-time. Under the 
assumption that the scalar field depends only on the coordinate y, the Einstein equations are 

-3aa" - 9{a'f + 9//^ 



a? 



a2 



where ' denotes the derivation with respect to y. The equation of motion of the scalar field is 
automatically satisfied by the solutions of because of Bianchi identity. 

For H > 0, the second equation of (0) gives an inequality 



J\2 „„// \ Tj2 



aa" >H'>0. (4) 

This inequality constrains the possible shape of a domain wall solution. At the peak of a 
domain wall defined by a' = 0, this inequality imposes a" < 0. Thus there can exist at most 
one peak, and a is a monotonically decreasing function of the distance from the peak. If we 
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assumed that a was bounded from below, a' and a" would approach zero in the limit \y\ —>■ oo, 
and would violate the inequality. Thus a must cross a = at finite y's on the both sides 
of a domain wall. The behavior of a near the vanishing point has a physical significance. 
Without tuning the behavior of a, (§) shows that the sources of the energy-momentum V{(j)) 
and (f)' diverge at the vanishing point and the point becomes a naked singularity. To avoid 
this situation the behavior of a near the vanishing point is constrained to have the form 

aiy) = -Hy + 0{y'), (5) 

where we are considering a domain wall whose peak is located at a negative value of y and 
that a vanishing point is located at ?/ = 0. Note that this behavior of a is similar to that of a 
Rindler horizon. In a series of the solutions presented in the vanishing points are horizons 
and the extended space-times are obtained. In this paper we will obtain and study another 
analytic solution of a de Sitter domain wall which is regular in the sense that the vanishing 
points have the behavior (|^). 

As discussed in [Q, the second order differential equations (^) can be rewritten in the 
following BPS-like first order differential equations by introducing a "superpotential" W{(j)): 

^ ^ 3mV(0) 

7 90 

a 2 



and the scalar field potential is given by 

Although this form of the differential equations is very convenient in the study of some general 
aspects of solutions and will be used later, another form discussed in is more convenient 
for the purpose of obtaining explicit solutions. With the substitution a = exp(— y4/2), the 
Einstein equations become 

V = -3A'^ + ^A" + 9H''exp{A), 

= 3A" -QH^exp{A). (8) 



"'^ However see [Q, where it is argued that solutions with singularities may also be physically meaningful in 
view of AdS/CFT correspondence. 



2 



A solution of these differential equations is a well-defined function of y. By solving ?/ as a 
function of A locally, one may regard the solution as a function of A rather than y with some 
care on multi- value. Now let us define g{A) by 



,2 3 dg{A) 



2 dA ' 

and integrate the second equation of (||) after multiplying A' . Then we obtain 



(9) 



g{A) = A" -AH'ei^^{A). (10) 

Here the constant shift ambiguity of the definition (^ of g{A) has been used to eliminate 
the integration constant in ([T0|). The g{A) is related to the previous form of the differential 
equations by 

g{A) = W{<py. (11) 
Thus the set of differential equations we will solve is given by 



l'2 



3dgiA) 



^ 2 dA ' 

A'^ = g{A) + 4H\xp{A), 

V ^ (12) 



For each choice of g{A) we obtain a solution ||^. We take the choice 

g{A) = -^{a^ + a), (13) 

where a and j3 are real constants and a = exp(— yl/2) as above. The first two equations of 
(H) give 

dY = 
da 



V a^-aa^ + f3^, (14) 

where Y = ^y, and we have chosen the signatures of the square roots corresponding to the 
right-hand side (larger y) region of the domain wall where increases and a decreases. The 
all over factor of a determines the length unit and we may fix it by setting a = 1 at the peak 
of the domain wall defined by a' = 0. Hence we take a = (3"^ — 1. The solutions to (0) are 
expressed by elliptic functions. Explicitly, we obtain 

a{y) = -sn{Hy,i/3'^), 
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where we have fixed the integration constants so that a vanishing point a = is located at 
1/ = and the scalar field takes = at the peak of the domain wall. The elliptic functions 
are defined by the inverse of the elliptic integral 

sn-i(z,A;)= / '^^ (16) 

io V(l -x2)(l-Px2) ' 

and cn('U, A;) = (1 — sn^(-u, A;))-^/^, dn('U, A;) = (1 — A;^sn^('U, A;))^/^. By solving a(0) and substi- 
tuting this into the third equation of (|T^), the potential is given by 



y(<p) = — + — ^2 — - cos (17) 




This scalar potential has the form of that of an axion with an instanton correction. 

By taking certain limits, our solution agrees with some of the known solutions. For H = 
cjS — >■ 0, the solution represents a static domain wall with H = and W{(f)) = 2csin(0/\/6) 
in (H) and (^. This "superpotential" was recently derived for a hypermultiplet in A/" = 2 
gauged supergravity derived from a non-homogeneous quaternionic space |l^ . This interesting 
development suggests a possible concrete route for the inclusion of our model into string- or 



M-theory, in spite of the widely applicable no-go theorem of |11|. For /3 = 1, the solution 



agrees with one of the solutions presented in p|. This can be checked by changing to the 
conformal coordinate dz = dy/a. For (3 ^ oo, the scalar field freezes at = and the 
solution is just for a five- dimensional cosmological constant [|T^. An analytic continuation to 
< does not work because 0(y) is ill defined. 

It is interesting to study how a static domain wall responds to an additional energy- 
momentum source. This kind of analysis was done for an energy-momentum source localized 



on a thin wall, and the standard cosmological evolution was reproduced |13, n] in the Randall- 
Sundrum approach [^. We can give a qualitative discussion on this question for a thick 
domain wall by interpreting judiciously our solution. Let us consider the change of the constant 
part of the potential (0) from that of the static case with W{4>) = 2csin(0/y6): 

Vi<P) = -^ + 5A + ^cos 
where 6 A is a positive constant. Comparing with ([T7|), we have 
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2/52 2 ■ ^ ' 



Solving the above equations we obtain 



= f . (20) 



Because of the redshift factor, the additional energy-momentum generated by 6A is well con- 
centrated near the wall, and hence 6A may be regarded as a source of the cosmological constant 
on the wall. Hence the relation (^) may be interpreted as ~ y^, which seems to be a good 
sign for the standard cosmological evolution for a thick domain wall. 

In the above analysis of obtaining (pOD, the regularity constraint is essential. If we 
allowed a naked singularity to appear, the expansion rate H could not be related to 6 A. In 
fact we can construct a solution with H = for any 6 A. For 7 = 1 [H = 0), is a first 
order differential equation for W{(j)): 



Since a' = at the peak of the wall, the trajectory of the solution goes through {(f), W) = (0, 0) 
in the {(f), W) plane. The solution W{(f)) = 2csin(0/v^) for the static case 5A = is a fine- 
tuned solution, and one can see that a slight change of the potential V{(f)) will result in a 
divergence of the asymptotic behavior of W{(f)) in large as observed in For most cases 
including the above modified potential, the trajectory of the solution in large enters the 
region where the potential energy V{(f)) is negligible in (PT]). Thus 



iy(0)~exph/-0 (22) 




for large 0. This region is where the scalar kinetic energy dominates over the potential energy. 
From (p2D and (j^), we have 

a{y) ~ (23) 

and this is a naked curvature singularity. Thus in our solution (|^), the curvature singularities 
are regularized by a de-Sitter expansion. This mechanism to use a de-Sitter expansion to turn 
a curvature singularity of a static soliton to a horizon has also appeared in the context of global 
U{1) vortex solutions [0 and in a codimension two non-supersymmetric soliton solution in 
IIB string theory |T7|. It would be very interesting to extend the analysis of the response of a 
thick domain wall to a more general energy-momentum source under the regularity constraint 
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The extension to the other cases than (0) would also be interesting. To avoid a curvature 
singularity at a = 0, the energy- momentum sources V and 0'^ should not diverge in the limit 
a —>■ 0. This is equivalent to the condition that both g{A) and dg{A)/dA do not diverge in 
the limit A —>■ oo. Thus a general choice g = CiO? + + C3 is also allowed for a regular 
solution, and we will obtain an explicit expression in terms of elliptic functions. In this case 
a non-trivial issue is whether the scalar field potential has an analytic expression in terms of 
(j) or not. 

Finally we will study the linear perturbations around our solution and prove its stability. 
The problem of obtaining the mass spectrum boils down to solving a Schrodinger equation 



+ Vt,e{z)]v = m'ip, (24) 



where Vt and K are the potential energies for the tensor perturbation [18| , 

ds'^ = a?{{j^i, + h^y)dx^dx'' + dz^), 

= V'' = 0, (25) 

and the scalar perturbation, 

ds^ = a^{{l + i)i)-i^^dx^'dx'' + (1 + i)2)dz^). (26) 
respectively [|T^, As for the tensor perturbation, the potential is given by 

V. ^ If + iH', (27) 

where 

(28) 

a dz 

The Schrodinger equation (|2^) can be rewritten in the form |^ 

gtQyp = mV, (29) 

where 

Q-i-\n. (30) 

Thus there is a normalizable state = c?!"^ with w? = 0, and the other states have > 0. As 
discussed generally in our solution is stable against the tensor perturbation. Substituting 
our solution ([T5|) into (^Tf), we obtain 

^' = W ^ " " ^""'^ ■ ^^^^ 
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This potential has a well near the domain wall peak and the low-energy states trapped in the 
well play the role as the graviton on the brane world. 



The potential for the scalar perturbation has the expression |18 



where {€) denotes the i-ih derivative with respect to z. Using (^), this potential can be 
expressed by a, W and 7. Assuming a form X = dia'jW + d2{a/'y)d'^W + d^a'y' /■y and after a 
straightforward tedious computation, we arrive at another expression 



K = + ^ + - 4^ ' (33) 

dz 6 



where X turns out to be 

X= \ ^ ' . 34 

dz 



The expression ( pBf ) shows that the Schrodinger equation for the scalar perturbation does not 
seem to have a non-negative expression like (|29|) . For if = it is evident that the remainder 
is non-negative, but for if > we cannot expect this in general. Especially, it is not so for 
our solution. Hence the question of stability whether there exist tachyonic states or not seems 
to need a case- by-case study. Substituting ( PTSD into the expression we obtain 

V; = 1^ (3/?2 + 5(1 - f?)Q' + a^) . (35) 



In the example studied in [0, the potential Vg is positive in the whole range of z and takes 
the maximum value at the domain wall peak. Thus it is evident that there does not exist 
any tachyonic or m = normalizable states in their example. In our example, however, the 
situation is more complicated. For > 3, there is a potential well with negative values near 
the domain wall peak. Hence it seems hard to answer the question of stability only from the 
shape of the potential. But fortunately, we have a useful relation 

^e = ^* + ^^- (36) 

Since we know that —d'^/dz^ + is non-negative, —d^/dz^ + Ve is positive definite. Thus our 
solution is stable against the scalar perturbation. 
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